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Abstract 

The non-Markovian features of three typical anomalous diffusing systems are studied by analyti- 
cally solving the generalized Langevin equation directly driven by three kind of internal structured- 
noises: harmonic noise, harmonic velocity noise and harmonic acceleration noise, respectively. The 
time-dependent reaction rate and the transmission coefficient are calculated by using of the reac- 
tive flux method. A startling behavior of Kramers-type effective reactive flow is witnessed in the 
harmonic noise and harmonic acceleration noise systems. 
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I. INTRODUCTION 



Substantial progresses have been witnessed during the last few years in understanding the 
phenomenon of anomalous diffusion in disordered media and the statistical thermodynamic 
cases that dose not obey the law of large number For example, the ergodicity breaking 

behavior accompanied by a nonunique stationary probability density for the corresponding 
embedded Markovian dynamics is of great interest both theoretically and experimentally 
4|-[7j]. This very feature is believed to originate from a vanishing effective friction of the 
system 8]. However, noticing that the friction is in general consanguineously related to 
the stochastic force (or noise) via the fluctuation-dissipation theorem (FDT) [9|, thus it is 
reasonable to argue that whether there exists some relationship between nonergodicity and 
the noise? Or in other words, is there a kind of noise environment that is relatively more 
favorable to the diffusing of a reactive particle? Such open problems are of great value to 
be considered. 

Recently, three types of thermal colored noises are proposed which provide three new 



typical structured-noise environments namely harmonic noise 
(HVN) and harmonic acceleration noise (HAN) respectively 
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N) , harmonic velocity noise 



llj . Obtaining from solving 



a second-order stochastic differential equation driven by a Gaussian white noise and so on, 
these noises are found to have very different types of power spectrums. The long time 
behavior of them have been carefully investigated showing some newfangled properties but 
how these noises affect the instantaneous dynamical process of a diffusing particle has not 
yet been intrinsically elucidated. The primary purpose of this paper is then to report our 
recent study on these problems focusing in particular on the instantaneous diffusing process 
of a particle in the structured-noise environments where a Kramers-type effective reactive 
flow is witnessed. 

The paper is organized as follows. In Sec. [Til the probability of passing over the barrier for 
a particle diffusing in the structured-noise environments is obtained by analytically solving 
the generalized Langevin equation (GLE). In Sec. Ill Ij. the time-dependent reaction rate is 
calculated by using of the reactive flux method and compared with the classical Kramers 
rate which is obtained in the case of Gaussian white noise. Sec. IIVI serves as a summary of 
our conclusion in which some implicit applications of this work are also discussed. 
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II. BARRIER PASSING PROBABILITY 



Structured-noise environments are generall y b elieved to originate from relevant kind of 



velocity-dependent system-reservoir couplings 



10Hl3|. Among the three types of structured 



noises, the simplest one is the harmonic noise (HN) which is a zero-mean Gaussian quasi- 
monochromatic noise having a power spectrum with a narrow Lorentzian peak centered at 
a finite frequency instead of zero. It is the solution of a second-order stochastic differential 
equation driven by a Gaussian white noise and has widely been used in many field of stud- 



ies 



14j-|l6|. The harmonic velocity noise (HVN), which can also be regarded as the time 



derivative of HN variable 



llj , can be obtained as a vortex transport of GLE in the presence 



of a mag netic field or a particle interacting via dipole coupling with a black-body radiation 
field [17| as well as in the presence of the velocity-coordinate coupling between the system 
and reservoir HQ . The harmonic acceleration noise (HAN) due to the velocity-velocity 



coupling in physical situations 
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13] cannot be procured exactly. However it can be taken 



as the second-order derivative of HN variable and the GLE driven by this noise can be solved 
numerically. 

For the physical systems which might exhibit an internal structured noise, the GLE 
describing the motion of a particle can be derived by using of the Zwanzig-Mori projection 



method hj, 



191 ] or the system-plus-reservoir method 20M22I] as 



mx{t) + / dt'/3(t - t')x(t') + d x U(x) = f (t), 



(1) 



where m is the mass of the particle, U(x) the potential energy, /3(t) the friction kernel due 
to the structured noise relating via Kubo's second kind of fluctuation-dissipation theorem 
(FDT II) 23] (£(£)£(£')) = mkBT(3(t — t') to the correlation function of each noise 
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(2b) 

(£(*)£(0>HVN, (2C) 



where r\ denotes the intensity of a Gaussian white noise, V and Q are the damping and 
frequency parameters in a RLC electric circuit respectively, fii and \ii are the roots of the 
characteristic equation /i 2 + T{i + Q 2 = 0. ks is the Boltzmann constant, T the temperature 
of the reservoir. 
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Due to the initially decoupling property between the system and reservoir, the structured 
noises are also of vanishing mean and not correlated with the initial velocity. The reduced 
distribution function of the particle can still be written as the classical Gaussian form 



W(x,t; x ,v ) 



1 



-exp 



21 



(3) 



V 2ira x (t) 

in which the average position (x(t)) and variance c^(t) can be obtained by Laplace solving 
the GLE. In the case of an inverse harmonic potential U(x) = — ^mu 2 x 2 , they are 



(*(*)) 



H(t')dt' 



x. 



dt x H(t - ti 



H(t)v 



(4) 
(5) 



where H(t) namely the response function can be yielded from inverse Laplace transforming 
H(s) = [s 2 + sj3(s) 
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with residue theorem [24| by inserting in the Laplace transfor- 



mation of the friction kernel function 



r?n 2 (s+r) 
r(s 2 + s r+Q 2 ) • 

sr]T 

s 2 + s r+n 2 ' 
sr](sr+n 2 ) 
r(s 2 +sr+n 2 ) • 



HN; 

HVN; 

HAN; 



(6) 



approaching from Eq.([21). 

Under the influence of thermal fluctuations, the probability density function (PDF) that 
is originally at one side of the barrier spreads with time, its center moves as well driven 
by its initial velocity The probability of passing over the saddle point (namely also the 
characteristic function) can then be determined mathematically by integrating Eq. 03) over 
x from zero to infinity 



X(x ,v ;t) 



W(x, t] Xq, V )dx, 



-erfc 
2 



(x(t)) 



(7) 



V2a x {t) 

Generally this will in the long time limit result in a finite real number range from to 1 
depending on the initial position and velocity of the particle with 1 for reactive trajecto- 
ries and for nonreactive ones. However, an fantastic x( x o,Va;t) is found in the case of 
structured-noise environments. 

As is shown in FIG. HJ where dimensionless parameters such as k^T = 1.0 are used in 
the calculation, the time-dependent barrier passing probability oscillates almost periodically 




FIG. 1: Time-dependent barrier passing probability x( x o, v o'it) f° r various kind of structured- 
noise environments: HN (solid), HVN (dash) and HAN (dot) plotted in the case of dimensionless 
parameters: m = ksT = V = rj = = 1.0 and VL 2 = 2.0 as well as xq = 0.0 and vq = 2.0. 

with time and that the "oscillating frequency" of which in the HAN case is about five times 
that of the HN or HVN case. This is a very fantastic phenomenon for the structured- 
noise environments that has never been witnessed before. The origin of this phenomenon 
may be retrospected to the implicit periodicity of the friction kernel function shown in Eq. 
02]). On this point one may argue: what a dynamical process on earth is the particle just 
experiencing; has a diffusing particle really passed the saddle point at t = t\ time? Many of 
such questions are now becoming wide open to dispute. 

However we noticed that although there lives some time- dependent uncertainties the 
periodical average of the barrier passing probability is near 0.5 for all the cases. This 
implies there would be on average about one-half amount of all the particles or a fifty-fifty 
chance for a particle to pass over the saddle point. Therefore the occurrence of an effective 
barrier escaping behavior is still expected for the particles diffusing in a structured-noise 
controlled metastable state potential. 
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III. RATE OF REACTIVE FLOW 



Now let us consider the time-dependent barrier passage of a particle diffusing in the 
structured-noise environments. Instead of solving the related Fokker- Planck equation [2? 
we derive the rate constant by using the method of reactive flux 
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271 ] . In the spirit of 



its calculations, the escape rate of a particle is defined by assuming the initial conditions to 
be at the top of the barrier, by investigating the ensemble of trajectories which start with 
identical initial conditions but experience different stochastic histories to determine whether 
a particle has successfully escaped or not. Mathematically it can be yielded from 



CO 



771 

k(t) = — - I dx I v W 8 t(xQ,Vo)S(x - x b )x(x ,v ;t)dv 



Qh 

k TST_^_ I v ^[-JpL\ x ( Xo = XblVo;t )do 0t (8) 



00 J —00 

771 f°° I TTlVn 



where Q is the partition function for reactions integrating over the distribution of the 
ground state; W s t(xo,v ) = H(xq,Vq) /ksT is an initial Boltzmann form stationary prob- 
ability distribution with H(xo,Vq) = \mvl + \7nuj\x\ representing the system Hamiltonian. 
k TST = -Qf^ e ~ VB ^ kllT is the classical rate resulted from the transition state theory (TST) 



28M30| in which Vb > denotes the height of the reactive barrier of the potential energy 
surface (PES). 

In order to make an isolated inspection on the dynamical corrections of k(t) to the 
TST rate, it is convenient to define the time-dependent transmission coefficient K(t) which 
describes the probability for a particle already escaped from the metastable well to recross 
the barrier point. By substituting Eqs.(H|, (JSJ and (GO) into Eq.([H]) we obtain 

r, ( 9) 



k B TE\t), 

which leads immediately to the Kramers formula for the rate constant derived from using a 



"flux over population" algorithm [3ll . |32j. However, a even funny behavior of nit) is found 
in the case of structured-noise environments. 

As is shown in FIG. [2] where the instantaneous variation of n(t) for various kind of 
structured-noise environments is plotted, we meet again in the long time with a periodical 
or more accurately semi-periodical oscillating n(t) in the HN and HAN case while a decaying 
to zero behavior in the HVN case. The stationary average of n(t) is about 0.43 in the HN 
case and 0.85 in the HAN case. This is also a startling phenomenon that has never been 




FIG. 2: Instantaneous variation of k(£) for various kind of structured-noise environments: HN 
(solid), HVN (dash) and HAN (dot). Identical dimensionless parameters are used as those shown 
in FIGfU 

witnessed before. This reveals that there exists on average an effective reactive flow in both 
HN and HAN noise environments. However, by comparing with those shown in Fig. [T] we 
can found that although a average barrier escaping probability may exist there is actually 
no effective reactive flow for a particle diffusing in the HVN noise environment. In other 
words, all the particles are absorbed in a finite range near the PES saddle point keeping 
on crossing and recrossing the barrier without contributing to an effective reactive flow. 
Therefore one may raise from these results that although some very features are believed to 
originate from the fluctuation and dissipation of the system, actually not all kinds of noise 
environments are beneficial to the diffusing particle, pessimistically some noise environments 
are even harmful. 



IV. SUMMARY AND DISCUSSION 

In summary we have studied in this paper the instantaneous dynamical process of a par- 
ticle diffusing in various structured-noise environments. The time- dependent barrier passing 
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probability, reactive rate and transmission coefficient are calculated in the framework of 
statistical Langevin reactive dynamics. A periodically oscillating barrier passing probability 
is witnessed in all the cases while an effective reactive flow only in the HN and HAN case. 
These meaningful results bring us with many totally brand-new understandings about the 
structured-noise environments. 

In further, due to their particular origins, this study will in no doubt renew our knowledge 
on some crucial problems in many field of studies such as the origin of nonergodicity in 
anomalous dissipation systems, reaction thermodynamics in chemical events, black-body 
radiation field studies and even those relevant to nanoscience experiments where the small 
systems are velocity-dependently coupled to a structured heat bath. 
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